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Abstract 

Highly specific spatiotemporal interactions between cognate molecular partners essentially sus- 
tain all biochemical transactions in the living matter. That such an exquisite level of accuracy 
may result from encountering forces solely driven by thermal diffusive processes is unlikely. Here 
we propose a yet unexplored strategy to experimentally tackle the long-standing question of a pos- 
sibly active recruitment at a distance of cognate partners of biomolecular reactions via the action 
of resonant electrodynamic interactions. We considered two simplified models for a preliminary 
feasibility investigation of the devised methodology. By taking advantage of advanced experimental 
techniques nowadays available, we propose to measure the characteristic encounter time scales of 
dually-interacting biopartners and to compare them with theoretical predictions worked out both 
in the presence or absence of putative long-range electromagnetic forces. 
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I. INTRODUCTION 



Living matter hosts a huge number of molecular players (i.e. proteins, nucleic acids) in- 
volved in simultaneous yet specific chemical reactions, despite an apparent lack of systematic 
spatial order. A phenomenological description of these biomolecular machineries at work 
often makes use of the concept of "recruitment", leaving usually unclear how biomolecules 
partners encounter or move toward their specific targets and sites of action. On this point 
significative progress has been made about DNA-protein interaction at short distance. This 
followed the puzzling problem posed by the E.Coli lac repressor-operator protein which was 
found to locate its specific DNA-binding site several orders of magnitude faster than the 
upper limit estimated for a diffusion-controlled process [l, 2]. A widely accepted approach 
to tackle this problem is the so-called facilitated diffusion, on which a vast literature exists 
(see for instance Refs.js-fj] and citations therein). To the contrary, for DNA-protein inter- 
actions and, more generally, for any dually-interacting biomolecules the mutual approach 
from a long distance is addresses to as 3D bulk diffusion and is not further studied (by " long 
distance" it is meant: much larger than the Debye screening length). Actually, at first in- 
spection, the mutual approach of cognate partners might well be driven by Brownian motion 
only, as at living temperature the ubiquitously distributed water molecules move chaotically 
in space, colliding with larger/heavier fluid components. On the latter, the total outcome 
of many simultaneous hits are forces of both random intensity and direction. Hence, by 
displacing themselves in a diffusive way through the inner cellular space, large molecules 
sooner or later will encounter their targets. 

A complementary proposal, which remains hitherto largely unexplored, is the possibility 
for molecules to interact at a distance via the electromagnetic field which is known to have 
sizeable magnitudes in living matter jl, 8]. In particular, electromagnetic attractive forces 
acting on a long range might, in specific conditions, facilitate the encounters of cognate part- 
ners, so that specific biomolecular reactions would occur more effectively than if dependent 
on stochastic motion only. Exploring this possibility, it should be stressed that the static 
dielectric constant e s of water is particularly high, e s ~ 80, at physiological temperatures. 
In addition to this dielectric screening, freely moving ions in the cellular medium tend to 
make the environment electrically neutral; accordingly, the Debye length in a biological 
environment is found to be smaller than ~ 10A, as was estimated on the basis of typical 
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ionic strength of the cytosol [6, 9J. Electrostatic interactions between electrically charged 
molecules at a distance larger than the Debye length are very unlikely. Conversely this is 
not necessarily the case for electro dynamic interactions [h], [ll| since the dielectric constant 
depends on the frequency of the electric waves under consideration. Among the latter, the 
interactions occurring between oscillating electric dipoles are of a particular interest since in 
many cases the long range nature of the interaction potential is essentially "activated" by 
the proximity of the dipole frequencies (resonance). In other words, two molecules whose 
dipole moments oscillate at the same frequency may undergo a so-called resonant interaction 
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131 ]. which is described by the potential U{r) oc — 1/r 3 with r the intermolecular distance 



(see Appendix). On the contrary, an off- resonance situation would produce a standard van 
der Waals-like potential, i.e. U(r) oc —1/r 6 , typically a short range interaction (see Ap- 
pendix). Such a frequency-selective interaction, when applied to a biological context, might 
be of utmost relevance during the approach of a molecule toward its specific cognate part- 
ner (s). To the best of our knowledge this proposition dates back to Jordan who advanced 
the idea that resonant interactions within a quantum framework could play a significant 
role in autocatalytic reactions or influence the process of biological synthesis in such a way 
that replicas of molecules present in the cell are formed [ijj]. His theory was questioned 
by Pauling pL5i], who estimated that such forces, supposed to occur only between identical 
molecules, could not be large enough to cause a specific attraction between proteins under 
the thermal conditions of excitation and perturbation prevailing in living organisms. Other 



attempts to explain 



der Waals forces 
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Diological selectivity have been made later on the basis of usual van 



171 ] . In parallel, in 1968, H. Frohlich proposed a dynamical model 18] 



to account for the capacity of biological systems to self-regulate, emphasizing that, under 
specific conditions of energy supply to these systems, part of this supply would not be totally 
thermalized but would be used to create order in response to environmental perturbations 
[19] ]. In particular, the normal polarization modes of a macromolecule (or of a part of it) 
may undergo a condensation phenomenon, characterized by the emerging of the mode of 
lowest frequency containing nearly all the energy supply [181 ] . Then, relying on this model, 



Frohlich suggested 
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J-]22j that - when occurring between two biomolecules - such dipole 
oscillations could be excited enough to overcome thermal noise leading to the above men- 
tioned frequency-dependent forces. Frohlich's seminal work has stimulated many theoretical 
investigations until our present days (see for example Refs. js, 23-25|). Moreover, a vast 
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literature is available about the experimental observation of low-frequency modes in the 



Raman and far infrared (TeraHertz) spectra of proteins 26( and DNA 27J. These spec- 
tral features are attributed to collective oscillation modes of the whole molecule (protein or 
DNA) or of a substantial fraction of its atoms. A-priori these collective oscillations of the 
molecular electric dipole moment could activate the mentioned long-distance attractive and 
selective recruitment interactions. However, a clear-cut experimental confirmation of the 
existence of the latter ones within a biological context at the molecular level is still lacking. 

In the present paper we consider a yet unexplored strategy to experimentally test, at 
least in simplified systems, whether these long-range recruitment forces are actually at work 
between typical actors of the broad variety of biomolecular reactions in living matter. On 
the basis of theoretical computations resorting to elementary and standard methods in the 
theory of stochastic processes on the one side, and recent progress on experimental methods 
on the other side, we make a first step toward the design of experiments to test whether 
such forces are actually at work in living matter. 

In Section [III we use two dynamical models to highlight qualitative and quantitative 
changes between Brownian and non-Brownian encounters of the macromolecular partners of 
a generic biochemical reaction. In Section IHIt we apply our models to the case of attractive 
electro dynamic potential U(r) oc — 1/r 3 expected to have effects at long distance, and then 
we report the numerical results that have been obtained with realistic parameters. Finally, in 
Section HVl we discuss how our findings can be used to design an experiment and we conclude 
that Fluorescence Cross Correlation Spectroscopy (FCCS) is an appropriate experimental 
tool to perform real-time measurements of the association kinetics of dually-interacting 
biopartners. It thus seems experimentally feasible to answer the basic questions formulated 
above by comparing the outcomes of the prospected experiments versus the theoretically 
predicted curves at different concentrations of the reactants. 

II. FIRST PASSAGE TIME MODELS 
A. Generalities 

Our idea is in principle a natural one: different kinds of forces must have different dynam- 
ical effects. Thus we attempted to devise an experimental protocol in order to discriminate 
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between the dynamics of purely random encounters between reaction partners versus en- 
counters driven by both a stochastic force plus a deterministic long-range force. Then, by 
experiments resorting on available techniques, we wondered whether it could be possible to 
discriminate between these different dynamical regimes. 

A natural way to proceed from the theoretical standpoint, that may be closely related 
to experimental as well as physiological conditions, is to consider an aqueous environment, 
initially containing N A particles of a species A and N B particles of a species B. Each 
molecule A is expected to interact with each molecule B in two ways : 

• As soon as the distance between A and B diminishes below a threshold 5, a biochemical 
reaction instantaneously takes place, so that the two molecules are not functional 
anymore and are considered as out of the system. 

• The particle A and the particle B interact at a distance via a two-body potential U(r) 
of an electrodynamic type, as long as the two molecules do not get closer than the 
distance 5. 

From a general point of view, the equations describing the dynamics of the system include 
both random and deterministic forces, and therefore may be given in the form: 

N B 

J^V A U(\v A , - r B>j \) + y/2i A kT$ A>i (t) 

3=1 

Na (1) 

^VBU(\r Ai -T BJ \) + y/2^itfZ Btj (t), 

i=l 

i = l,...,N A , and j = 1, . . . , N B . 
Here, m A , m B correspond to the masses, r A ^, r B j to the positions, and 7^4, j B to the 
friction coefficients of the constituents of each species. T stands for the temperature in 
the solution, and k is the Boltzmann constant. £(£) is the random process modeling the 
fluctuating force due to the collisions with water molecules, usually represented as a Gaussian 
white noise process for which (£ Ai {t)^ k {t')) = 5 a p5 ik 8(t-t'), where a, = 1,2,3 are related 
to each component of the £^j's. The same relation is valid also for the $, B /s- 

A-priori, equations describe a very complex dynamics, even in the absence of ran- 
domness. For example, assuming that the potential is, for each pair of molecules, of the 



m A - 



d 2 r A , 
dt 2 



m B - 



d 2 r Bt 
dt 2 
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dr 



A A 



dt 



-7s- 



dt 
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form U(r) = c\jr m — C2/r n , with C\,C2 constants, m,n G N, m > n and n < d (long-range 
condition if d is the spatial dimension), the Hamiltonian subset of this system is actually a 
nonlinear classical iV-body system whose phase space is entirely filled with chaotic trajecto- 



ries 
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At this stage, the addition of random forces may imply that the representative point 
of the system nontrivially wanders in phase space, despite the presence of dissipative terms 
which, in principle, would generate trivial attractors. Indeed, in the overdamped limit, when 
the acceleration terms can be neglected, one is dealing with a randomly perturbed first-order 
nonlinear dynamical system which, as integrability is exceptional, is expected to display a 
complex (chaotic) dynamics. Nevertheless, instead of undertaking the numerical integration 
of Eqs.Q, we decided to look, as a first step, for some analytic result that can be obtained 
at the cost of some simplification of the system. 

Because the reaction between two particles A and B occurs the first time they come 
sufficiently close together, we will have to focus on first passage times of a simplified version 
of system (0Q). Generally, first passage or first return time statistics are difficult to examine 
in dynamical systems, then leading one to model the system under study by keeping only 



29 



Here we 



its salient characteristics, either in a deterministic or stochastic manner 
rather choose to still work with equations (1), but to reduce drastically the dimensionality 
of the system then leading to keep in the model under study its salient characteristics only. 
This was achieved by noting that Eqs. (JT|) describe the mutual interaction between the two 
sets of particles, A and B, but neither the A nor the B particles interact among themselves. 
The trajectories of the A particles are indirectly coupled only through the dynamics of the B 
particles. Thus, as a first simplifying hypothesis, we assumed that the B particles are fixed, 
and as a consequence the dynamical behaviors of the different A particles are independent. 

This leads to the decoupling of the individual equations in ([1]) and hence to the introduc- 
tion of a one-dimensional model representative of the generic dynamics of a single A particle. 
This model is considered below according to two different versions and solved according to 
standard methods [31]. 



B. Model 1: absorbing plus reflecting boundaries 

Let us consider one fixed molecule B located at the position z = and one molecule A, 
initially located at z = x (see Figure 1). We first suppose that if A reaches the boundary 



6 



z = L of the domain, it is reflected back to z < L; whereas when A reaches the position 
z = 5 for the first time, it is absorbed. The random trajectory z(t) of the molecule A may 
be given, as previously, in the form 



dz 

-jl = v » 

dv (2) 
m^ = - 1 v + F(z) + £(t). 

For times much larger than the characteristic time m/7, equations (J2]) will then relax to a 
state in which dv/dt — > 0. This approximation is justified by the fact that the biomolecules 
involved in reactions of interest (protein-protein or DNA-protein) typically weigh thousands 
of Daltons, and thus the characteristic relaxation times in aqueous medium are very short. 
Therefore, equations ((2]) for the A particle can be simplified as 



! = £W + fe (t) . (3 ) 

at 7 V 7 



3l| 



As it is well known, the one dimensional Langevin initial value problem 

dz 

— = a(z,t) +b(z,t)£(t) , z(t ) = x, 

is equivalent to the Fokker-Planck equation (FPE) for a probability p(z,t\x,to) of finding 
the particle at z at time t, given it was at x at to ^ t 

d d Id 2 

0jjp(z> t\ x , to) = - ^[ a ( z > *M Z > t\ x , *o)] + 2 H z ' f ) 2 P( z ^ t\ x > f o)] ■ ( 4 ) 
From ()3]), one thus obtains 

^P(z, t\x, t ) = ~-]^ [ F ( z )p( z i t\x, t )} + k ^ q z2 p(z, t\x, t ), (5) 
which is also known as the Smoluchovski equation. 

We now look at the time T at which the reaction between A and B occurs. That is, 
the first time when particle A reaches z — 8. Since we are considering an absorbing barrier 
at z = 5 and a reflecting barrier at z — L, the probability P(T > t) and the one that the 
particle would still be in the interval [8, L] at time t are the same 
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FIG. 1: (Color online) A generic initial condition of Model 1 (t = 0). Here x is the initial distance 
between A and B; 5 is the distance at which A and B react, and L is the position of the reflecting 
barrier for A and the position of B is fixed. 

P(T>t) = ^ dzp{z,t\x,0) := G(x,t). 

Besides, as F(z) and kT do not explicitly depend on t, p(z, t\x,0), and thus G(x,t), are 
homogeneous processes, such that 

G(x,t) = f dzp(z,0\x,-t). (6) 

J 8 

This implies that G(x, t) satisfies the same partial differential equation oip(z, 0|x, —t) for 
z fixed, that is, a backward Fokker-Planck equation 



9 , n , s l n/ N 9 , nl , fcT9 2 

— p(z, 0|x, -t) = -F(x)— p(2, O x, -t) H — 

at 7 ax 7 ax- 



leading to 

d ltd d 2 1 

-G(x,i) = — |F(x)— G(x,i) - kT—G(x,t)\. (7) 

Here, the initial condition p(z, 0\x, 0) = 5(x — 2) (here 5 is the Dirac functional) clearly 
gives 



G(x, 0) = 1, if 5 < x < L; and G(x, 0) = 0, if not, (8) 

whereas the absorbing condition at 5 and the reflecting boundary condition at L allow 
to write, respectively 



d 

G(S,t) = and T^G(x,t) 
ox 



= ,Vf>0. (9) 

x=L 



If one focuses on the mean first passage time r(x), which represents a characteristic time 
scale of the reaction, one has by definition 



Tlx) 



OO 

t-F(T<t)dt 



d 

t— G(x,t)dt = I G{x,t)dt, 



(10) 



after integration by parts. Then, by integrating Eq.flT]) between t = and t = oo, and 
using the fact that G(x, 0) = 1 and G(x, oo) = 0, we find that t(x) must satisfy the following 
ordinary differential equation 



1 1 [F(x) dT{x) hT d2r{x) 
dx dx 2 
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with boundary conditions t(S) = dr(x)/dx\ x= L = 0, as it follows from equations ([9]) and 
(|T0|) . The solution is found to be [3l| 



T(x) = l dy W)l dz ^ {z) ' 



with 



ip(x) = exp 



F(s) 



ds 



since F(x) 



is kT 

-dU(x) / dx. This gives for r{x): 



exp 



U(x) - U(S) \ 
kT J 



T{x) = kfJ s dyGXP \kT J 



kT 



(11) 



(12) 



It can easily be checked that the mean first-passage time in presence of an attracting 
deterministic potential, generically written as U(x) oc — x~ n with a given n > 0, is smaller 
than the mean first-passage time with Brownian motion only, i.e., when U = 0. Since 
exp (x~ n ) is a decreasing function of x, we can find an upper limit for the second integral 
and thus get 

-L 



Tlx) < 



J_ 

kT 



dy exp 



1 



kTy n J J y 



dz exp 



kTy r 



J_ 
kT 



dy / dz := r(x 
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More explicitly 



r(x) Bwn = ^f[(L- 5? ~ (L - xf] = ^(x - 6)(2L - 6 - x). (13) 
C. Model 2: two absorbing boundaries 

Let us now consider the alternative model where two particles B are fixed at positions 
z = and z = I, so that the particle A, initially located at z = x (see Figure EJ), is absorbed 
as soon as it reaches z = 5 or z = I — 5. Such a model is mathematically similar to the 
previous one. 
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FIG. 2: (Color online) A generic initial condition of Model 2 (t = 0). Two molecules B are fixed 
at the boundaries x = and x = I; x and Z — x are the initial distances between A and the Bs and 
5 is the distance at which A and B react. 

Simply the deterministic force and the boundary conditions have to be modified. Conse- 
quently Eq.(([7|)) is to be replaced by 



^G(x, t) = ~[ l F (x) ~ HI - *)] J^(z, t) - kT^G(x, t) | (14) 

i-l-S 

now G(x, t) is defined as G(x, t) — / <iz t|x, 0). The initial condition p(z, 0|x, 0) 
<5(x — z) gives 



G(ac,0) = l, if 5 and G(x,0)=0, if not, (15) 

and the absorbing boundary conditions give 

G(6,t)=0 and G(l-8,t) = ,Vt>0. (16) 

The mean first-passage time r(x) defined above, then satisfies the ordinary differential 
equation 
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-l = i([F(a;)-F(/-x)]^-A;T rf2r(x) 



7 [ dx dx 2 

with boundary conditions 



T (S) = T (l -6) = 0. 



The solution is found to be 



3l| 



t(x) = 2 



U 



Z-<5 I f x 1 f w J 



with 



iix) =exp >l ns)- F(l-s) d 1 = ^ f £,(*) - - x) - U(S) + - *) 



since _F(a;) = —dU(x)/dx. After simplification, one has 



;i7) 



where 



(sj = exp 



U(s)-U(l-s) 



kT 

Similarly to the Model 1, the expression for the mean reaction time with Brownian motion 
only (U = 0), is particularly simple 

T (x) B ™ = ^L{x-S){l-6-z). (19) 

To summarize, in this Section we have obtained the general form of the mean first-passage 
time t(x), that is the average time needed by molecule A to reach the molecule B (or one 
molecule B in the case of Model 2), as a function of the initial intermolecular distance x and 
temperature T, for both Model 1 [Eq. (lT2]) ] and Model 2 [Eq. (II8]) ]. respectively. The same 
function is given for randomly driven encounters between the reaction partners in Eq. (|T3|) 
for Model 1, and in Eq.(JTSJ for Model 2. 
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III. QUANTITATIVE THEORETICAL PREDICTIONS 



In order to answer the question of whether it would be feasible to experimentally detect 
the possible existence of a deterministic attractive force through which the cognate partners 
of biochemical reactions interact at long distance, we first have to delimit the physical 
context, choose the domain of physical parameters and provide the analytic form of the 
two-body interaction potential. In what follows, a long range resonant potential potential 
U(x) = — C/x 3 is considered. As discussed in the Introduction, this kind of interaction 
can have sizeable effects at long distances at variance with London - Van der Waals 1/x 6 
interactions (see also the Appendix). Following Frohlich 22], a lower bound for the coefficient 
of this potential is given by C ~ he 2 (Z A Z B ) 1/2 /2Mu e'(u ), where Z A and Z B denote the 
number of charges of averaged mass M and charge e contributing to the dipole moment of the 
molecules A and B respectively, whereas ujq stands for their oscillation frequency; e'(uo) is the 
real part of the dielectric constant of the interposed medium. In particular, it is interesting 
to remark that in the expected range of oscillation frequencies for the setup of collective 
dipole oscillations in macromolecules (that Frohlich estimated to be ujq ~ 10 11 — 10 12 Hz) the 
value of £'(^0) drops down to a few units 



32J thus allowing a much smaller screening of the 



interactions with respect to the static case. In this context, we use Za = Z B = 1000 (see 
Ref. |8|) and the proton mass for M. A convenient unit system remains to be chosen. For 
the numerical tabulations of t(x), instead of c.g.s. units we use /xm, kDa, and /xs, with the 
following definitions 1 /xm =10~ 4 cm, 1 kDa = 10~ 21 g and 1 /xs = 10 _6 s. 

In this system of units we evaluate the lower bound of C which is found to be C ~ 
10 _30 erg.cm 3 = 0.1 kDa./xm 5 ./xs~ 2 . Henceforth, we shall consider C varying from 0.1 to 
10 kDa./xm 5 ./xs -2 . These values are given with a degree of arbitrariness that can be reduced 
by considering that C = 10 corresponds to the physical situation where U(x) ~ kT at 
x ~ 0.1/xm. Hence the choice C G [0.1, 10] is a very cautious estimate with respect to those 
existing in the literature about a possibly larger range of action (it has been surmised by 
rohlich and others that U(x) might become comparable with kT at x ~ 1/xm or more 
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81, 1231, 



241]). Among other constants appearing in equations (1121) and (1181) . the friction 



coefficient 7 of the molecule A has been estimated according to Stokes' law 7 = Q^rj{T)lZ, 
where rj(T) corresponds to the viscosity of water at temperature T and 1Z stands for the 
hydrodynamic radius of the molecule. The value of 1Z has been set equal to 5 • 10~ 3 /xm which 
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is the typical diameter of a biomolecule with a mass in the interval 50 — lOOkDa (proteins 
and DNA fragments); the same value has been fixed for the reaction radius 5 introduced in 
both models : 1Z = 8 — 5 • lCT 3 /xm. 

All the computations of r(x) have been performed by means of MATLAB programs. Also, 
as MATLAB does not allow to perform direct integrations over non-rectangular domains, 
integrals with variable limits in equations (fl2|) and (fl8|) have been first "vectorized" [33J] for 
each x to calculate t(x) with a recursive adaptive Simpson quadrature (MATLAB quadl 
function). Further checks on the reliability of the method have been done through direct 
numerical integration of the Langevin equation by means of a standard Euler-Heun algo- 
rithm and by averaging over 10 4 different realizations of the random walk. Minor precision 
problems especially when x ~ 5 have been thus detected and corrected in what follows. 



A. Model 1 

We have computed r(x) and r(x) Bwn by means of Eqs. (fT2]) and (IT31 respectively, where 
we have set R = 5 = 5 • lCT 3 /zm and U (x) = — C/x 3 , as detailed above. The position L of the 
reflecting barrier characteristic of Model 1 has been fixed so as x <C L for all x. In particular, 
L = 10/zm and a maximal value for x equal to 1/im have been used. Figure [3] displays the 
numerically found shapes of both functions r(x) and T(x) Bwn computed at T = 300K and 
for different values of the attractive potential coefficient C. A first check on the reliability 
of the plotted results is done by observing that t(x) < t(x) Bwti for all the x values while 
both curves merge at large x values, as expected when the resonant attraction is wiped out 
by thermal noise. In particular, as we always considered i < i, the asymptotic behavior 
of r(x) Bwn is then proportional to x as required by Eq. fflBl . On the contrary, at smaller 
x, t(x) Bwu bends downwards to slightly smaller values with respect to the extrapolated 
linear dependence (this happens when x is no longer much larger than 5). At variance, 
the pattern of t(x) has two asymptotic limiting behaviors: at large x values it joins the 
Brownian curve r(x) Bwn , and at small x values it is r(x) ~ x 5 (a power-law characteristic 
of the 1/x 3 form of the potential); the latter might be anticipated on the basis of simple 
dimensional arguments since the l.h.s. of Eq.([3]) has the dimensions [dz/dt] = It' 1 while the 
r.h.s. leads to [C/x 4 ] = [C]l~ 4 in a purely deterministic regime. By combining the two, for a 



generic time scale [r] = t associated with a displacement length x we get r ~ [x] [361 ] . The 
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two limiting behaviors are bridged by a steep transition pattern which moves rightward or 
leftward according to the value of C, as shown on Figure [3] the stronger the potential the 
larger the x-values at which t(x) displays the knee joining the x 5 functional dependence. 
The transition pattern of r(x) is steep since the reflecting barrier is located far from the 
only one molecule B. 

In any case, r(x) Bwn is found to exceed r(x) by a factor of, say, 10 at definitely larger 
a;- values and, what is more relevant, at longer values of first encounter time : with C = 0.1, 
such a difference occurs at x ~ 300A where r Bwn ~ 6ms and r ~ 600/is, while with C = 1.0, 
it occurs at x ~ 640A where r Bwn ~ 10ms and r ~ 1ms, and with C = 10 at x ~ 1400A 
where r Bwn ~ 20ms and r ~ 2ms. As we shall see in the next Section, should we interpret x 
as the average distance between any two reacting molecules in three dimensions, this range 
of x- values (between a few hundreds Angstroms and 1/im) is easily attained by varying 
the concentrations of the reactants between a few micro-Moles down to one nano-Mole. 
Notably, the encounter times belong to an interval of values easily accessible by means of 
optical detection methods. 

A priori, further qualitative indications on the possible presence of attractive deterministic 
forces between cognate partners could be observed by modifying the temperature of the 
system. As shown in Figure HI t(x) and r(x) Bwn plotted for three different values of T 
confirm that the x 5 functional dependence is purely deterministic as T has no influence 
within this domain. On the contrary, r(x) Bwn displays the same dependence on T for all 
values of x . Surprisingly, the steep transition pattern of r(x) at intermediate values of x is 
characterized by a temperature dependence which is inverted compared to the Brownian case 
: in presence of an intermolecular potential, the higher the temperature the larger the first- 
passage time of A at x = 5. Finally note that the temperature range considered in Figure 
H] is a broad one (T = 200, 300, and 400K). Nevertheless, as the physiological temperature 
range corresponds only to a few percent around 300K, it is likely that variations of the first 
passage time at different temperatures are too weak to be experimentally detectable within 
such an interval. In particular, computations performed for temperature differences of 10K 
(typically 290, 300, 310K) with the Model 1 show variations of r less than five percent of its 
value in the Brownian well as in the case of Brownian plus deterministic force. 
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FIG. 3: (Color online) Model 1: mean encounter time r between two molecules A and B, initially 
placed at a distance x one from the other. Dotted lines are asymptotic behaviors. Dashed line 
refers to purely random encounters. Solid line refers to the combined effect of a random force plus 
a deterministic one derived from the potential U(x) = — C/x 3 . Figure (a) refers to C = 0.1. Figure 
(b) refers to C = 10. 
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FIG. 4: (Color) Model 1: temperature dependence of t(x) for C = 1.0. Red solid and dashed 
curves refer to T = 200K, green dashed and solid curves refer to T = 300K, blue dashed and solid 
curves refer to T = 400K. 

B. Model 2 

We have also plotted r(x) and r(x) Bwn computed according to Eqs. (lT8]) and ( Tl~9l) respec- 
tively, where R — 6 — 5 • 10 _3 /im, I = 2x, and U(x) = —C/x 3 as in the case of Model 1. 
Figure O displays the numerically found shapes of both functions r(x) and r(x) Bwn computed 
again at T = 300K and for different values of the attractive potential coefficient C. 

At first check, the characteristics of r(x) prevailing for Model 1 (i.e., the x 5 functional 
dependence at small x, and the tendency of r(x) to join r(x) Bwn at large x) also apply 
here. Nevertheless, the asymptotic x-dependence of r(x) Bwn is now proportional to x 2 as 
obtained from equation Eq.(fl9l) by replacing I by 2x when x ^> 5. On the other hand, the 
absence of a reflecting barrier makes the steep transition feature for t(x) disappear and be 
replaced by a mild crossover at intermediate values of x. This steeper pattern for the first 
passage time as a function of x - in the case of purely Brownian diffusion - entails a less 
pronounced separation between r(x) and r(x) Bwn . In fact, a separation between these two 
curves by a factor of 10, (to make the same kind of comparison that we did for Model 1) 
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FIG. 5: (Color online) Model 2: mean encounter time t between two molecules A and B, initially 
placed at a distance x from each other. Dotted lines are asymptotic behaviors. Dashed line refers 
to purely random encounters. Solid line refers to the combined effect of a random force plus a 
deterministic one derived from the potential U(x) = —C/x 3 . Figure (a) refers to C = 0.1. Figure 
(b) refers to C = 10.. 
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FIG. 6: (Color) Model 2: temperature dependence of t{x) for C = 0.1. Red solid and dashed 
curves refer to T = 200K, green dashed and solid curves refer to T = 300K, blue dashed and solid 
curves refer to T = 400K. 

occurs for C = 0.1 at x ~ 200A for which r Bwn ~ 1.1/is and r ~ 0.1/xs; while for C = 1.0 
this happens at x ~ 600A for which r Bwn ~ 30/xs and r ~ 3/xs; and, finally, for C = 10 at 
x ~ 11 70 A for which r 5 ^™ ~ 120/xs and r ~ 12/is. Therefore, we can see that Model 2 is 
more constraining than Model 1, in the sense that, at equal values of C (that is, at equal 
strength of the long range interaction) smaller intermolecular distances and a much faster 
tracking of the dynamics of the reactants are needed to discriminate with the same degree 
of confidence (arbitrarily set as a factor of 10) between random and non random encounters 
of the reactant molecules. 

The temperature dependence of both r(x) Bwn and t(x) is reported in Figure El The 
main features of t(x) in Model 1 are likewise present in Model 2. In addition, we can see 
that the inversion in the temperature dependence with respect to the Brownian case (which 
was characteristic of the steep transition pattern of Model 1) is no longer there. Thus, the 
peculiar temperature dependence of this steep transition pattern could be mostly attributed 
to the presence of the reflecting barrier characteristic of Model 1. Likewise in Model 1, we 
used again T = 200, 300 and 400K, even though computations carried out at physiological 
temperature again yield too weak variations of t(x) to be experimentally detectable. 
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IV. DISCUSSION AND CONCLUDING REMARKS 



The numerical results reported in the preceding Section are in favor of a positive answer 
to the main question addressed by the present work. In fact, the numerical study of Models 
1 and 2 revealed qualitative differences in the mean first passage time r between the case 
of a pure Brownian diffusion of the molecule A (see Section [IT]) and the case in which an 
attractive (resonant) potential U is added to a random force. In particular, in the latter 
case, the functional dependence of r on the initial distance x between the molecule A and 
its target (molecule(s) B) demonstrates the existence of different patterns in the two models 
depending on the range of the x values considered : 

• a deterministic pattern at small x values (small initial separations), characterized by 
a power law representative of the potential under consideration (x 5 for the resonant 
potential used, x p+2 for a general potential of the form U(x) oc x~ p ); 

• a Brownian pattern at large initial separations, proportional to x or to x 2 depending 
upon the symmetry of the system (x for an asymmetric situation as described by 
Model 1, and x 2 for a symmetric one as described by Model 2); 

• a steep transition pattern joining the two asymptotic ones in the case of Model 1, and 
a smooth crossover joining the two asymptotic ones in the case of Model 2. 

Although complementary computations revealed some interesting features in the tem- 
perature dependence of r, the corresponding degree of variation in a laboratory accessible 
interval of temperature is too weak to be experimentally detectable. 

In any case, it is obvious that x must constitute an experimentally accessible control 
parameter so that the results mentioned above may be used to predict the possible role 
of long-range intermolecular forces in biological processes. Notably, such an approach is 
not so usual. Indeed, most of the attempts made hitherto in this direction have resulted in 
experimental measurements of association constants k a (characteristic of a reaction medium), 
which are predictable from the Smoluchowski theory also when intermolecular forces are 
considered 
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35| . The focus of Smoluchovski theory is on the association constant k a 



which represents the probability for two molecules to react per time unit, irrespective of their 



position. In the case of Brownian encounters, this is given by 



4 fcf = AixRD = 4n5kT/~f 
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where R is the reaction radius that can be approximated to 5 from the current study, and 
D is the sum of the diffusion coefficients of the two cognate partners. In the presence of 
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some interaction potential U, one has k* = 4nR*D where R has been replaced by 
R* — R I j r - 2 e - u i r )/ kT d r j . Now, if k a is experimentally measured for some reaction and 

it turns out that k^ < k a , then this would indicate that some deterministic force is in action 
but one can hardly find out the law of the interaction potential because after integration 
over r there is no one-to-one correspondence between R* (thus k*) and the functional form 
of U(r) [24|. On the other hand, in measuring k a < k^ , we cannot be sure that the reaction 
is simply diffusion-driven because, in this case, chemical times could be long enough to 
make k a smaller than the corresponding Brownian value. The advantage of our dynamical 
approach is that our models still apply by choosing 5 as the distance at which A and B get 
in contact without reacting, and with the experimental technique discussed below (FCCS) 
we can make a distinction between the association time and the chemical times. 

In the present situation x values might to some extent be considered in three dimension 
as the average distance between two molecular partners A and B, while this quantity can be 
easily controlled in laboratory experiments by varying the concentrations of the reactants. 
Given the concentrations Ca = N a /Va and Cb = N b /Vb (with Va,b = the initial volumes 
and Na,b = the number of molecules of the two species respectively; remark that these num- 
bers are controlled through the molarity, i.e. a definite fraction of the Avogadro number), 
we get the estimate x = C~^ 3 for the average intermolecular distance from the average 
concentration C av = (Na + Nb)/V, where the reaction volume V = Va + Vb- In practice, 
as an example, with C av = InM we have x ~ 1/im as the average distance between any two 
molecules, while with C av = 1/iM we have x ~ 1000A. By working at equimolarity, that is 
Ca = Cb, then C~^ 3 is a good estimate of the average distance between one A and one B 
molecule. Working with nano-Moles of DNA and proteins (enzymes, transcription factors) 
is quite standard in molecular biology experiments. With such concentrations of reactants, 
both models (1 and 2) predict that the first passage time - that can be interpreted as the 
average encounter time between one A and one B molecule initially located at intermolecular 
distances of a few thousands of Angstroms - varies in the interval between a few tens of mi- 
croseconds to about one millisecond in the presence of an attractive deterministic force that 
would sum up to the random force. On the contrary, in the very same conditions, random 
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only driven encounters would exceed the above mentioned encounter times by one or two 
orders of magnitude. Again, the distance at which sizeable differences could be observed 
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ue of the resonant potential parameter 



24| suggest that these long-range reso- 



may vary significantly depending on the actua 
C. On the other hand, estimates in literature js, 
nance interactions could be effective up to distances in the order of 1/im (the action range 
is estimated by computing the distance at which the resonance interaction energy equals 
the level of thermal noise kT). In this respect, in the preceding Section we have limited 
ourselves to cautious estimates for the parameter C, focusing on conservative assumptions 
for average encounter time varying in the interval ICT 5 — I CP 3 seconds which can be readily 
detected with the aid of Fluorescence Cross-Correlation Spectroscopy (FCCS technology). 
This is a powerful technique which is being increasingly applied to the study of diffusion 
and chemical reaction rates in complex biological systems using fluorescently labeled macro- 
molecules 
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39(. FCCS measures the spontaneous fluctuations of fluorescences 8Fi(t) and 



<5i<2(t) that arise from the diffusion of fluorescently labeled molecules of type 1 and 2, re- 
spectively - illuminated by two laser light beams of different colors - into or out of an open 
sampling volume. Even though the size of the detection volume is diffraction limited, the 
autocorrelation functions of 5F 12 {t) and the cross-correlation function (5Fi(t)5F 2 (t)) can be 
altered by processes occurring on smaller spatial scales. These correlation functions provide 
information on diffusion properties of fluorescent molecules. 

Of course, we are well aware of the fact that the models studied here are simplified de- 
scriptions of the reality. Indeed, protein-protein and protein-nucleic acid interactions in vivo 
generally take place within complex structural scaffolds such as the membrane cytoskeleton 
or the chromatin envelope, which are themselves the subject of highly dynamical regula- 



tions (e.g., Refs. |40|, |4l|); and may a 



of the given reactions (e.g., Refs. 
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so possibly interfere with the spatiotemporal control 



43]). Should resonant electrodynamic interactions be 



involved within such an intricate context, it seems illusory at this stage to assess realistic C 
values simply based on the proposed experiments. In fact, regarding protein-DNA in vivo 
(physiological) interactions for instance, it may well be that the putative values fluctuate de- 
pending on a host of variables, possibly including - in a non-mutually exclusive way, charges 
on proteins and DNA, the effect of surrounding electrolytes, the nucleic-/amino-acid com- 
positions, the length of accessible DNA, etc. However, we stress here that our initial goal, as 
described in this article, is to merely probe whether or not biological partners can take ad- 
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vantage, besides thermic diffusion, also of long-distance (0.1-1/zm) forces of electrodynamic 
origin to eventually interact. If established, this novel concept would then in turn open new 
avenues of research to investigate long-standing biological issues, e.g., on the precise defi- 
nition of which variables exactly pertain on protein-DNA interactions, and how a diffusing 
protein particle may actually recognize the particular cognate DNA site among many other 
locations also available. Since we do not expect dramatic qualitative changes out of the 
numerical simulations of Eqs. ([1]) in three dimensions [44j , an experimental setup providing 
a practical realisation of what has been investigated in the present work could be devised 
by resorting - as experimental probes - to three broad classes of interactions: protein-DNA, 
protein- RNA, and protein-protein (ligand- receptor). As DNA and RNA molecules have 
not a preassigned length, it is implicitly understood that only short fragments are to be 
considered (some tens or a few hundreds of base pairs, that is, oligonucleotides or plas- 
mides respectively). The proteins interacting with DNA or RNA can be processing enzymes 
(helicases, polymerases, recombinases) or transcription factors normally bound at promot- 
ers, enhancers, insulators, or silencers. Thus, for example, one could choose two molecular 
species consisting, respectively, of a short double stranded DNA molecule (for example a 
synthetic oligonucteotide of ~ 100 base pairs or even less) and a protein with a site specific 
affinity for the chosen DNA molecule (i.e., a transcription factor). By labeling the DNA 
molecules and proteins with standard fluorophores their dynamical behavior can be followed 
by means of FCCS microscopy at different concentrations C = Ca — Cb of the reactants to 
get a characteristic time scale as a function of x = C^ 1 ^ 3 . In this way such an experimental 
set up should provide - after data fitting - an estimate of the constant C for the resonant 
potential considered above. Thus, C = would mean that the reactants meet only under the 
action of Brownian diffusion, whereas C / would prove the existence at the same time of 
the long-range interactions evoked throughout this paper and give quantitative information 
about them. 
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V. APPENDIX 



For the sake of clarity and to help the reader to get a hold of the physical origin of the 
U(R) oc — 1/R 3 potential referred to throughout the present work, this Appendix provides 
some theoretical elements about the interaction of oscillating electric dipoles. 

To begin with, let us recall some basic fact on this subject. Two atoms (or two small 
molecules) A and B in their ground states with no net charge excess and vanishing average 
dipole moment (i.e. both are nonpolar) interact through the London - Van der Waals 
dispersive force. The origin of this interaction is as follows. Though the expectation values 
of the dipole operators are zero for nonpolar atoms, quantum fluctuations are responsible 
for their instantaneous non vanishing dipolar moments. This entails a non zero dispersion of 
dipole moment operator. The energy of the two isolated atoms is corrected at first order by 
the dipole-dipole interaction energy which is proportional to the average dipole moments, 
thus it vanishes when both atoms are in their ground states. Instead, the second order 
perturbative correction, due to the coupling between instantaneous dipole fluctuations, is 
found to be proportional to 1/R 6 . (In a QED framework the London - Van der Waals 
interaction stems from the exchange of virtual photons between the atoms). This is a short 
range potential, so called because the exponent of the power law of R is strictly larger than 
3, the dimension of physical space. London - Van der Waals interactions are of generically 
weak intensity, whereas they likely become of prime importance in a biological context when 
acting at short distances (below the Debye length) together with additional interactions of 
chemical type 45]. 

Remarkably, the first order perturbative correction may be non- vanishing under a degen- 
eracy condition. Indeed, if one or both atoms are in an excited state, provided that the 
condition for exchange symmetry is fulfilled, that is, they have common eigen-energies in 
their spectra, it can be shown 46( that the interaction energy is now proportional to 1/R 3 , 
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a long range potential. 

Interactions of similar kinds to those just mentioned between two atoms (or small 
molecules) could exist between macromolecules with an oscillating electric dipole moment. 
In this case, the oscillating dipole moment would not be due to the electron motions but, 
rather, to conformational vibrations. As already mentioned in the Introduction, this in- 
teraction between the oscillating electric dipole moments of reacting macromolecules could 
play a relevant role in living matter. In fact, as already quickly recalled in the Introduction, 
the high static dielectric constant of water together with the considerable amount of ions 
present in living cells tend to screen any electrostatic interaction beyond a distance of a few 
Angstroms. However, this electrostatic opaqueness does not hold for an oscillating field: the 
higher the frequency of an oscillating field the more transparent an aqueous salted medium. 



In fact, the value of the die 
function of the frequency 



ectric constant of water at room temperature is a decreasing 
32] and, for example, already at ITHz e(oj = 10 12 ) ~ 4; likewise, 
the imaginary (dissipative) part of the dielectric constant (which is proportional to the con- 
ductivity of the medium due to the presence of free ions) is inversely proportional to the 



frequency of the oscillating electric field (according to the Drude equation [48|), so that at 
suitably high frequency can be negligible. 

Let us a now study the basic mechanism of interaction between two oscillating electric 



471 ] . As we show below, these 



dipoles before discussing its application to biomolecules 
oscillating dipoles can activate long-range forces that will be shown to be frequency selective. 
We consider a one dimensional simplified model in which the dipoles oscillate at frequencies 
ua and ub respectively. Then a computation of the interaction energy between A and B 
can be given which, despite the simplified treatment, allows to grasp some basic physical 
facts. 

Let fx a and lib be the masses of the two oscillators, let their dipole moments be parallel 
and given by qZxrA and qZsTB, and assume that their mutual separation R is such that 
R ^> tai^b-i then we can write the interaction Hamiltonian as 

H = + & + ^ w + ^ M + TB ' m 

where Zi, i — A, B, stands for an effective number of charges which account for the average 
value of the dipole moment of the oscillator i; ( is a geometrical factor depending on the 
orientation of the dipoles with respect to the line joining them (on which the distance R is 
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measured). Then, introducing a mean mass M defined so that ha = MZa and /ib — MZb, 
the Hamiltonian becomes 



H= ^M ^ + P ^ + \ Mu) ^ + \ Muj B r B + Jp r A r B , 



(21) 



where the transformations (Zi) x l 2 ri — > r\ and pi/iZi) 1 ^ 2 — >■ p^, % — A,B, have been 
introduced (the variables rj and pj are still canonically conjugated) and we put (3 = 
(q 2 (ZAZ B ) 1 ^ 2 /47re , where e is the dielectric constant of vacuum - in the absence of a 
material medium between the oscillators - to be replaced by e{oj) when a medium is present. 
In matrix form this also reads 

/ co 2 (3 /MR 3 \ /~.\ 



K (3 /MR 3 



UJ, 



I 



ta 



(22) 



c 



Matrix C is real and symmetric, thus diagonalizable by means of an orthogonal transfor- 
mation. Let oj 2 + and u 2 __ the eigenvalues of C (homogeneous to squared frequencies). Under 
the action of this transformation the Hamiltonian can be cast in the form of the sum of two 
decoupled oscillators, that is, 



1 11 

H = (p 2 + p 2 ) + -Mu 2 ,r 2 + -Mtu 2 r 2 , 

2M v + ; 2 + + 2 



(23) 



and it can be easily shown that 



V2 



K 2 + - B 2 )±{K 2 -c B 2 ) 2 + J^} 



1/2' 



1/2 



(24) 



By considering ta, tb, Pa and ps as observables subject to standard commutation relations, 
the energy values of the system are obviously given by 



E hu: + ( n + + - j + fvjj- in + - 



(25) 



where n + , n_ e N, that is, are integers. Now, let us consider two opposite physical situations 
depending on the relative values of the frequencies oua and w B of the oscillators. 



1. Consider oja ^b (or, equivalently, oja «C wb), we have 



V2 



(u A 2 + Ub) ± K 2 - ^b) \ 1 + 



4/3 



{u 2 - u 2 ) 2 M 2 R e 



1/2 



1/2 



(26) 
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and the denominator of the last term is large enough to give at the lowest order 
expansion 

1/2 



1 

7! 



K 2 + Ui) ± K 2 - Ui) ± K 2 _ 2 ^ 2)M2fi6 + 



2/3 



1/2 



■ P 2 



where ua,b stands for uja in the computation of u + and ujb in the computation of w_. 
By substituting this expression in Eq. fT25l) we get 

(27) 



£ = ftu^ ( n + + - J + ftw B ( n_ + - I + 



2(w A 2 - u 2 )M 2 R & \u A V 2/ w B V ~ 2, 
The first two terms correspond to the unperturbed energies of the oscillators A and 
B considered as isolated (R — > oo) while the last term provides the lowest order 
correction to the unperturbed energy of the system and due to the interaction, this 
interaction potential energy is proportional to R~ 6 . Note that this is functionally the 
same as the London - Van der Waals interaction but of a remarkably different physical 
origin (real oscillations instead of quantum fluctuations). 

2. To the contrary, at resonance, that is, ua — ^b — the eigen-frequencies 021]) are 
simply given by 

-W^siW- (28) 

At long distances (imposed by the reality condition for u± in this equation) we can 
develop cu± near uj and replace such a development into Eq. (l2"5"j) to obtain 

hd hd 2 

E = fiw(n + + „_ + 1) + (n+ - „_) - (n + - „_ + 1) + ... (29) 

The first order correction to the energy of the system corresponds to the interaction 
energy between the two oscillators at resonance and is proportional to R~ 3 . If both 
oscillators are in their ground states, i.e. n + = n_ = 0, the first contribution to 
the interaction energy in Eq. (l29|) vanishes as well as the force given above. The first 
non vanishing term is again proportional to i?" 6 . But if the lowest of these modes 
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(oj-) gets more excited than the other (oj + ) then the consequence is the activation 
of an attractive long-range frequency-selective force. A repulsive force could also be 
activated in case n + > n_. 



In the context of Frohlich's theory [8|, |18|, |49|, |50[ the above described mechanism of 



resonant interaction between oscillating dipoles was surmised to have a great poten- 
tial relevance for fundamental biological processes at the molecular level. Frohlich 
proposed a model describing the coupling between the elastic vibrations of macro- 
molecules and the resulting time variations of their dipole moment; the model predicts 
that one or a few Fourier modes of the dipole field oscillation should be strongly (co- 
herently) excited provided that the energy supply rate exerted on the macromolecule 
by its environment exceeds a threshold value. This energy supply is assumed to de- 
pend on the biological activity of the environment (metabolic energy). The strongly 
excited mode of oscillation of the molecular dipole moment should be due a collective 
oscillation either of the entire molecule or of a subgroup of its atoms. The consequence 
of such collective oscillations would be to activate selective long-range recognition and 
attraction between cognate macromolecular partners via the above described mech- 
anism of resonant interaction. Experimental evidence of the existence of collective 
excitations in macromolecules of biological relevance is available for polynucleotides 
(DNA and RNA) j^] and for proteins 26] in the Raman and far infrared (TeraHertz) 
spectroscopic domains. 
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